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1 The RG: what, why and how 

Imagine that you have some problem in the form of a partition function 

Z{a,b) ^ f dx f dye-'^("'+^')e-''("+«)' (1) 



where a, b are the parameters. 

First consider & = 0, the gaussian model. Suppose that you are just inter- 
ested in X, say in its fluctuations. Then you have the option of integrating 
out y and working with the new partition function 



Z{a) = N j dxe-"'' (2) 

where N comes from doing the y-integration. We will ignore such an x- 
independent pre-factor here and elsewhere since it will cancel in any averaging 
process. 

Consider now the nongaussian case with b ^ Q. Here we have 



Z{a',b'...) = / dx 



dxe-"' '■^\-b'x^-c'xO+... 

where a' , b' etc., define the parameters of the effective field theory for x. These 
parameters will reproduce exactly the same averages for x as the original ones. 
This evolution of parameters with the elimination of uninteresting degrees of 
freedom, is what we mean these days by renormalization, and as such has 
nothing to do with infinities; you just saw it happen in a problem with just 
two variables. 

The parameters &, c etc., are called couplings and the monomials they 
multiply are called interactions. The x^ term is called the kinetic or free-field 
term. 
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Notice that to get the effective theory we need to do a nongaussian in- 
tegral. This can only be done perturbatively. At the simplest tree Level, we 
simply drop y and find b' = b. At higher orders, we bring down the non- 
quadratic exponential and integrate in y term by term and generate effective 
interactions for x. This procedure can be represented by Feynman graphs in 
which variables in the loop are limited to the ones being eliminated. 

Why do we do this? Because certain tendencies of x are not so apparent 
when y is around, but surface to the top, as we zero in on x. For example, we 
are going to consider a problem in which x stands for low-energy variables 
and y for high energy variables. Upon integrating out high energy variables 
a numerically small coupling can grow in size (or initially impressive one 
diminish into oblivion), as we zoom in on the low energy sector. 

This notion can be made more precise as follows. Consider the gaussian 
model in which we have just a ^ 0. We have seen that this value does not 
change as y is eliminated since x and y do not talk to each other. This is 
called a fixed point of the RG. Now turn on new couplings or "interactions" 
(corresponding to higher powers of x, y etc.) with coefficients &, c and so 
on. Let a', 6' etc., be the new couplings after y is eliminated. The mere fact 
that b' > b does not mean b is more important for the physics of x. This is 
because a' could also be bigger than a. So we rescale x so that the kinetic 
part, x^, has the same coefficient as before. If the quartic term still has a 
bigger coefficient, (still called b'), we say it is a relevant interaction. E b' < b 
we say it is irrelevant. This is because in reality y stands for many variables, 
and as they are eliminated one by one, the coefficient of the quartic term will 
run to zero. If a coupling neither grows not shrinks it is called marginal. 

There is another excellent reason for using the RG, and that is to un- 
derstand the phenomenon of universality in critical phenomena. I must re- 
gretfully pass up the opportunity to explain this and refer you to Professor 
Michael Fisher's excellent lecture notes in this very same school many years 
ago p. 

We will now see how this method is applied to interacting fermions in 
d ~ 2. Later we will apply these methods to quantum dots. 



2 The problem of interacting fermions 

Consider a system of nonrelativistic spinless fermions in two space dimen- 
sions. The one particle hamiltonian is 

H^^-t^ (4) 

where the chemical potential ^ is introduced to make sure we have a finite 
density of particles in the ground state: all levels up the Fermi surface, a 
circle defined by 
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Kl/2m^fi (5) 

are now occupied and occupying these levels lowers the ground-state energy. 

Notice that this system has gapless excitations above the ground state. 
You can take an electron just below the Fermi surface and move it just above, 
and this costs as little energy as you please. Such a system will carry a dc 
current in response to a dc voltage. An important question one asks is if this 
will be true when interactions are turned on. For example the system could 
develop a gap and become an insulator. What really happens for the d = 2 
electron gas? 

We are going to answer this using the RG. Let us first learn how to do 
RG for noninteracting fermions. To understand the low energy physics, we 
take a band of of width A on either side of the Fermi surface. This is the 
first great difference between this problem and the usual ones in relativistic 
field theory and statistical mechanics. Whereas in the latter examples low 
energy means small momentum, here it means small deviations from the 
Fermi surface. Whereas in these older problems we zero in on the origin in 
momentum space, here we zero in on a surface. The low energy region is 
shown in Figure ^ 




Fig. 1. The low energy region for nonrelativistic fermions lies within the annulus 
concentric with the Fermi circle. 



To apply our methods we need to cast the problem in the form of a path 
integral. Following any number of sources, say |2I we obtain the following 
expression for the partition function of free fermions: 

Zo = J diPdbpe^" (6) 

where 

So = j d^K p duoi^iuo, K) (^lu: - i^^!^-^^ V^(c^, K) (7) 

where ip and ip are called Grassmann variables. They are really weird objects 
one gets to love after some familiarity. Most readers can assume they are 
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ordinary integration variables. The dedicated reader can learn more from 
Ref. la. 

We now adapt this general expression to the annulus to obtain 

Zo = J di^dApe'^" (8) 

where 

Sq^ d9 duj dh^P{iu; - V k)^. (9) 

Jo J-oo J~A 

To get here we have had to approximate as follows: 

^\^K^.k^.,k (10) 

where k — K ^ Kp and vp is the fermi velocity, hereafter set equal to unity. 
Thus A can be viewed as a momentum or energy cut-off measured from the 
Fermi circle. We have also replaced KdK by Kpdk and absorbed Kp in ^ 
and ip. It will seen that neglecting k in relation to Kp is irrelevant in the 
technical sense. 

Let us now perform mode elimination and reduce the cut-off by a factor s. 
Since this is a gaussian integral, mode elimination just leads to a multiplica- 
tive constant we are not interested in. So the result is just the same action 
as above, but with |fc| < A/s. Let us now do make the following additional 
transformations : 



(w',fc') = s{uo,k) 

(^'(c', k')3\u\ k')) = s-^l\i,{-, -), V^(-, -)). (11) 

s s s s 

When we do this, the action and the phase space all return to their old 
values. So what? Recall that our plan is to evaluate the role of quartic inter- 
actions in low energy physics as we do mode elimination. Now what really 
matters is not the absolute size of the quartic term, but its size relative to 
the quadratic term. Keeping the quadratic term identical before and after 
the RG action makes the comparison easy: if the quartic coupling grows, it is 
relevant; if it decreases, it is irrelevant, and if it stays the same it is marginal. 
The system is clearly gapless if the quartic coupling is irrelevant. Even a 
marginal coupling implies no gap since any gap will grow under the various 
rescalings of the RG. 

Let us now turn on a generic four-Fermi interaction in path-integral form: 

54 = / ^(4)^^(3)^(2)^(1)7.(4, 3, 2, 1) (12) 



where / is a shorthand: 
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A 

dki I diOi (13) 



— oo 



At the tree level, we simply keep the modes within the new cut-off, rescale 
fields, frequencies and momenta , and read off the new coupling. We find 



u{k\uj' ,6) ^u{k' /s.uo' /s,e) (14) 

This is the evolution of the coupling function. To deal with coupling con- 
stants with which we are more familiar, we expand the functions in a Taylor 
series (schematic) 

U = Uo + kui + fc^M2... (15) 

where k stands for all the fc's and w's. An expansion of this kind is possible 
since couplings in the Lagrangian are nonsingular in a problem with short 
range interactions. If we now make such an expansion and compare coeffi- 
cients in Eqn. H14I) . we find that uq is marginal and the rest are irrelevant, 
as is any coupling of more than four fields. Now this is exactly what happens 
in 04, scalar field theory in four dimensions with a quartic interaction. The 
difference here is that we still have dependence on the angles on the Fermi 
surface: 

liO = m(6'i,6'2,6'3,6'4) 

Therefore in this theory we are going to get coupling functions and not a 
few coupling constants. 

Let us analyze this function. Momentum conservation should allow us to 
eliminate one angle. Actually it allows us more because of the fact that these 
momenta do not come form the entire plane, but a very thin annulus near 
Kp- Look at the left half of Figure El Assuming that the cutoff has been 
reduced to the thickness of the circle in the figure, it is clear that if two 
points 1 and 2 are chosen from it to represent the incoming lines in a four 
point coupling, the outgoing ones are forced to be equal to them (not in their 
sum, but individually) up to a permutation, which is irrelevant for spinless 
fermions. Thus we have in the end just one function of two angles, and by 
rotational invariance, their difference: 



m(0i, 02, 01, 02) = i^(0i - O2) = F{e). (16) 

About forty years ago Landau came to the very same conclusion^] that a 
Fermi system at low energies would be described by one function defined on 
the Fermi surface. He did this without the benefit of the RG and for that 
reason, some of the leaps were hard to understand. Later detailed diagram- 
matic calculations justified this picture The RG provides yet another way 
to understand it. It also tells us other things, as we will now see. 

The first thing is that the final angles are not slaved to the initial ones if 
the former are exactly opposite, as in the right half of Figure 2. In this case, 
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Fig. 2. Kinematical reasons why momenta are either conserved pairwise or re- 
stricted to the BCS channel. 

the final ones can be anything, as long as they are opposite to each other. 
This leads to one more set of marginal couplings in the BCS channel, called 



The next point is that since F and V are marginal at tree level, we have 
to go to one loop to see if they are still so. So we draw the usual diagrams 
shown in Figure 3. We eliminate an infinitesimal momentum slice of thickness 
dA8Xk = ±A. 

These diagrams are like the ones in any quartic field theory, but each 
one behaves differently from the others and its its traditional counterparts. 
Consider the first one (called ZS) for F. The external momenta have zero 
frequencies and lie of the Fermi surface since oj and k are irrelevant. The mo- 
mentum transfer is exactly zero. So the integrand has the following schematic 
form: 



The loop momentum K lies in one of the two shells being eliminated. Since 
there is no energy difference between the two propagators, the poles in lo lie 
in the same half-plane and we get zero, upon closing the contour in the other 
half-plane. In other words, this diagram can contribute if it is a particle-hole 



u(0i,-0i,(?3,-f?3) = vie^ - Oi) = v{e). 



(17) 




(18) 
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Fig. 3. One loop diagrams for the flow of F and V. The last at the bottom shows 
that a large momentum Q can be absorbed only at two particular initial angles 
(only one of which is shown) if the final state is to lie in the shell being eliminated. 

diagram, but given zero momentum transfer we cannot convert a hole at 
to a particle at +A. In the ZS' diagram, we have a large momentum transfer, 
called Q in the inset at the bottom. This is of order Kp and much bigger 
than the radial cut-off, a phenomenon unheard of in say cf)'^ theory, where all 
momenta and transfers are bounded by A. This in turn means that the loop 
momentum is not only restricted in the direction to a sliver dA, but also in 
the angular direction in order to be able to absorb this huge momentum Q 
and land up in the other shell being eliminated (see bottom of Fig. ©. So 
we have du ~ dt^, where dt = dA/A. The same goes for the BCS diagram. 
Thus F does not flow at one loop. 

Let us now turn to the renormalization of V. The first two diagrams are 
useless for the same reasons as before, but the last one is special. Since the 
total incoming momentum is zero, the loop momenta are equal and opposite 
and no matter what direction K has, —K is guaranteed to lie in the same shell 
being eliminated. However the loop frequencies are now equal and opposite 
so that the poles in the two propagators now lie in opposite half-planes. We 
now get a flow (dropping constants) 




(19) 
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Here is an example of a flow equation for a coupling function. However by 
expanding in terms of angular momentum eigenfunctions we get an infinite 
number of flow equations 



one for each coefficient. These equations tell us that if the potential in angular 
momentum channel m is repulsive, it will get renormalized down to zero ( a 
result derived many years ago by Anderson and Morel) while if it is attractive, 
it will run off, causing the BCS instability. This is the reason the y s are not 
a part of Landau theory, which assumes we have no phase transitions. This 
is also a nice illustration of what was stated earlier: one could begin with a 
large positive coupling, say V3 and a tiny negative coupling W5. After much 
renormalization, W3 would shrink to a tiny value and W5 would dominate. 

3 Large- iV approach to Fermi liquids 

Not only did Landau say we could describe Fermi liquids with an F function, 
he also managed to compute the response functions at small lo and q in terms 
of the F function even when it was large, say 10, in dimcnsionless units. Again 
the RG gives us one way to understand this. To this end we need to recall 
the the key ideas of "large-N" theories. 

These theories involve interactions between N species of objects. The 
largeness of N renders fluctuations (thermal or quantum) small, and enables 
one to make approximations which are not perturbative in the coupling con- 
stant, but are controlled by the additional small parameter 1/A^. 

As a specific example let us consider the Gross-Neveu modeljB] which is 
one of the simplest fermionic large- theories. This theory has N identical 
massless relativistic fermions interacting through a short-range interaction. 
The Lagrangian density is 



Note that the kinetic term conserves the internal index, as does the in- 
teraction term: any index that goes in comes out. You do not have to know 
much about the GN model to to follow this discussion, which is all about the 
internal indices. 

Figure 01 shows the first few diagrams in the expression for the scattering 
amplitude of particle of isospin index i and j in the Gross-Neveu theory. The 
"bare" vertex comes with a factor X/N. The one-loop diagrams all share a 
factor \^ /N'^ from the two vertices. The first one-loop diagram has a free 
internal summation over the index k that runs over N values, with the con- 
tribution being identical for each value of k. Thus, this one-loop diagram 




(21) 
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Fig. 4. Some diagrams from a large-N theory 



acquires a compensating factor of which makes its contribution of order 
A^/-/V, the same order in 1/N as the bare vertex. However, the other one- 
loop diagrams have no such free internal summation and their contribution 
is indeed of order 1/iV^. Therefore, to leading order in 1/N, one should keep 
only diagrams which have a free internal summation for every vertex, that 
is, iterates of the leading one-loop diagram, which are called bubble graphs. 
For later use remember that in the diagrams that survive (do not survive), 
the indices i and j of the incoming particles do not (do) enter the loops. 
Let us assume that the momentum integral up to the cutoff A for one bubble 
gives a factor —11 (A, qext), where qext is the external momentum or frequency 
transfer at which the scattering amplitude is evaluated. To leading order in 
large- the full expression for the scattering amplitude is 



r{qext) = 



1 



A 



(22) 



N 1 + \n{A, qext) 

Once one has the full expression for the scattering amplitude (to leading 
order in \/N) one can ask for the RG flow of the "bare" vertex as the cutoff 
is reduced by demanding that the physical scattering amplitude F remain 
insensitive to the cutoff. One then finds (with t = \n{Ao/A)) 



dF(qext) 



dX 



= -A- 



.dn{A,qext) 



(23) 



dt dt dt 

which is exactly the flow one would extract at one loop. Thus the one-loop RG 
flow is the exact answer to leading order in a large- iV theory. All higher-order 
corrections must therefore be subleading in 
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3.1 Large-N applied to Fermi liquids 

Consider now the ipip — tpif) correlation function (with vanishing values of 
external frequency and momentum transfer). Landau showed that it takes 
the form 

where _Fb is the angular average of F{9) and xo is the answer when F = Q. 
Note that the answer is not perturbative in F . 

Landau got this result by working with the ground-state energy as a 
functional of Fermi surface deformations. The RG provides us with not just 
the ground-state energy, but an effective hamiltonian (operator) for all of 
low-energy physics. This operator problem can be solved using large iV- 
techniques. 

The value of N here is of order Kp/A, and here is how it enters the 
formalism. Imagine dividing the annulus in Fig. into N patches of size 
{A) in the radial and angular directions. The momentum of each fermion 
is a sum of a "large" part {0{kp)) centered on a patch labelled by a patch 
index i = 1, ...N and a "small" momentum {0{A) within the patch|2]. 

Consider a four- fermion Green's function, as in Figure Q). The incoming 
momenta are labelled by the patch index (such as i) and the small momentum 
is not shown but implicit. We have seen that as yl 0, the two outgoing 
momenta are equal to the two incoming momenta up to a permutation. At 
small but finite A this means the patch labels are same before and after. 
Thus the patch index plays the role of a conserved isospin index as in the 
Gross-Neveu model. 

The electron-electron interaction terms, written in this notation, (with 
k integrals replaced by a sum over patch index and integration over small 
momenta) also come with a pre- factor of (~ A/Kp)- 

It can then be verified that in all Feynman diagrams of this cut-off the- 
ory the patch index plays the role of the conserved isospin index exactly as 
in a theory with N fermionic species. For example in Figure Q in the first 
diagram, the external indices i and j do not enter the diagram (small mo- 
mentum transfer only) and so the loop momentum is nearly same in both 
lines and integrated fully over the annulus, i.e., the patch index k runs over 
all N values. In the second diagram, the external label i enters the loop and 
there is a large momentum transfer {0{Kp)). In order for both momenta in 
the loop to be within the annulus, and to differ by this large q, the angle of 
the loop momentum is limited to a range 0{A/ Kp). (This just means that if 
one momentum is from patch i the other has to be from patch j. ) Similarly, 
in the last loop diagram, the angle of the loop momenta is restricted to one 
patch. In other words, the requirement that all loop momenta in this cut-off 
theory lie in the annulus singles out only diagrams that survive in the large 
N limit. 
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The sum of bubble diagrams, singled out by the usual large- consider- 
ations, reproduces Landau's Fermi liquid theory. For example in the case of 
X, one obtains a geometric series that sums to give x = i+Fo ' 

Since in the large N limit, the one- loop /3- function for the fermion-ferniion 
coupling is exact, it follows that the marginal nature of the Landau parame- 
ters F and the flow of V , Eqn. (|2()|l . are both exact as yl ^ 0. 

A long paper of mine Ref. ( 2 ) explains all this, as well as how it is to be 
generalized to anisotropic Fermi surfaces and Fermi surfaces with additional 
special features and consequently additional instabilities. Polchinski 'S' inde- 
pendently analyzed the isotropic Fermi liquid (though not in the same detail, 
since it was a just paradigm or toy model for an effective field theory for 
him). 



4 Quantum dots 

We will now apply some of these ideas, very successful in the bulk, to two- 
dimensional quantum dots[3inj- tiny spatial regions of size L ~ 100 — 200nm, 
to which electrons are restricted using gates. The dot can be connected weakly 
or strongly to leads. Since many experts on dots are contributing to this 
volume, I will be sparing in details and references. 

Let us get acquainted with some energy scales, starting with Z\, the mean 
single particle level spacing. The Thouless energy is defined as Et = ?i/T, 
where r is the time it takes to traverse the dot. If the dot is strongly coupled 
to leads, this is the uncertainty in the energy of an electron as it traverses 
the dot. Consequently the g (sharply defined) states of an isolated dot within 
Et will contribute to conductance and lead to a (dimensionless) conductance 
n — 

The dots in question have two features important to us. First, motion 
within the dot is ballistic: L^i, the elastic scattering length is the same as 
L, the dot size, so that Et = where t^iT- is the Fermi velocity. Next, 

the boundary of the dot is sufficiently irregular as to cause chaotic motion 
at the classical level. At the quantum level single-particle energy levels and 
wavefunctions (in any basis) within Et of the Fermi energy will resemble 
those of a random hamiltonian matrix and be described Random Matrix 
Theory (RMT) 8 . We will only invoke a few results from RMT and they 
will be explained in due course. 

At a generic value of gate voltage Vg the ground state has a definite num- 
ber of particles N and energy £jv. If £n+i — £n = aeVg {a is a geometry- 
dependent factor) the energies of the N and A^-l-l-particle states are degener- 
ate, and a tunnelling peak occurs at zero bias. Successive peaks are separated 
by the second difference of £n, called A2, the distribution of which is mea- 
sured. Also measured are statistics of peak- height distributions [TU| . [TTl [T^ . 
which depend on wavefunction statistics of RMT. 

To describe the data one needs to write down a suitable hamiltonian 
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Hu ^'^^a^'^ ^ 2 ^ Vapfsciclc^cs (25) 

(where the subscripts label the exact single particle states including spin) and 
try to extract its implications. Earlier theoretical investigations were confined 
to the noninteracting limit: V = and missed the fact that due to the small 
capacitance of the dot, adding an electron required some significant charging 
energy on top of the energy of order it takes to promote an electron by one 
level. Thus efforts have been made to include interactions |lHI IT^ 1151 llfij . 

The simplest model includes a constant charging energy UoN'^/2 jllil 17]. 
Conventionally Uq is subtracted away in plotting A2- This model predicts a 
bimodal distribution for A2: Adding an electron above a doubly-filled (spin- 
degenerate) level costs C/q -I- s, with s being the energy to the next single- 
particle level. Adding it to a singly occupied level costs Uq. While the second 
contribution gives a delta-function peak at after Uq has been subtracted, 
the first contribution is the distribution of nearest neighbor level separation 
e, of the order of A. But simulations [TS| and experiments^! E| produce 
distributions for A2 which do not show any bimodality, and are much broader. 

The next significant advance was the discovery of the Universal Hamilto- 
nian ^^El- Here one keeps only couplings of the form Vafd/3a on the grounds 
that only they have a non-zero ensemble average (over disorder realizations). 
This seems reasonable in the limit of large g since couplings with zero aver- 
age are typically of size 1/g according to RMT. The Universal Hamiltonian 
is thus 

Hu^Y: ^"4..c., + - +^(Y<A.i) (Ywa) (26) 

a,s ^ a '' ^ [3 

where s is single-particle spin and S is the total spin. The Cooper coupling A 
does not play a major role, but the inclusion of the exchange coupling J brings 
the theoretical predictions^! 1151 ^ into better accord with experiments, 
especially if one-body "scrambling" ^! 120! 1211 and finite temperature 
effects are taken into account. However, some discrepancies still remain in 
relation to numerical^B! and experimental results^! at > 2. 

We now see that the following dot-related questions naturally arise. Given 
that adding more refined interactions (culminating in the universal hamilto- 
nian) led to better descriptions of the dot, should one not seek a more system- 
atic way to to decide what interactions should be included from the outset? 
Does our past experience with clean systems and bulk systems tell us how to 
proceed? Once we have written down a comprehensive hamiltonian, is there 
a way to go beyond perturbation theory to unearth nonperturbative physics 
in the dot, including possible phases and transitions between them? What 
will be the experimental signatures of these novel phases and the transitions 
between them if indeed they do exist? These questions will now be addressed. 
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4.1 Interactions and Disorder: Exact results on the dot 

The first crucial step towards this goal was taken by Murthy and Mathur 
\'2l^ . Their ideas was as follows. 

- Step 1: Use the clean system RG described earlier (eliminating mo- 
mentum states on either side of the Fermi surface ) to eliminate all states 
far from the Fermi surface till one comes down to the Thouless band, that 
is, within Et of Ep. 

We have seen that this process inevitably leads to Landau's Fermi liquid 
interaction (spin has been suppressed): 



where 6*, 9' are the angles of k, k' on the Fermi circle, and Um is defined 



A few words before we proceed. First, some experts will point out that 
the interaction one gets from the RG allows for small momentum transfer, 
i.e., there should be an additional sum over a small values q in Eqn. H27|l 
allowing k — > k + q and k' ^ k' — q. It can be shown that in the large g 
limit this sum has just one term, at g = 0. Unlike in a clean system, there is 
no singular behavior associated with q — > and this assumption is a good 
one. Others have asked how one can introduce the Landau interaction 
that respects momentum conservation in a dot that does not conserve 
momentum or anything else except energy. To them I say this. Just think 
of a pair of molecules colliding in a room. As long as the collisions take 
place in a time scale smaller than the time between collisions with the 
walls, the interaction will be momentum conserving. That this is true for 
a collision in the dot for particles moving at vp, subject an interaction of 
range equal to the Thomas Fermi screening length (the typical range) is 
readily demonstrated. Like it or not, momentum is a special variable even 
in a chaotic but ballistic dot since it is tied to translation invariance, and 
that that is operative for realistic collisions within the dot. 
- Step 2: Switch to the exact basis states of the chaotic dot, writing the 
kinetic and interaction terms in this basis. Run the RG by eliminating 
exact energy eigenstates within Et- 

While this looks like a reasonable plan, it is not clear how it is going to be 
executed since knowledge of the exact eigenfunctions is needed to even write 
down the Landau interaction in the disordered basis: 




(27) 



by 




(28) 



Vcp^s = iE^iO- 9') C(k)<^Mk') - rai^'m^) 



kk' 



X [0^(k')05(k)-0^(k)05(k')] 



(29) 
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where k and k' take g possible values. These are chosen as follows. Consider 
the momentum states of energy within Ex of Ep. In a dot momentum is 
defined with an uncertainty Ak ~ 1/L in either direction. Thus one must 
form packets in k space obeying this condition. It can be easily shown that g 
of them will fit into this band. One way to pick such packets is to simply take 
plane waves of precise k and chop them off at the edges of the dot and nor- 
malize the remains. The g values of k can be chosen with an angular spacing 
2Tr/g. It can be readily verified that such states are very nearly orthogonal. 
The wavefunction (j)5(k) is the projection of exact dot eigenstate S on the 
state k as defined above. 

We will see that one can go a long way without detailed knowledge of the 
wavefunctions (f>s(k-)- 

First, one can take the view of the Universal Hamiltonian (UH) adherents 
and consider the ensemble average (enclosed in <> ) of the interactions. RMT 
tells us that to leading order in 1/g, 

<C(ki)</';(k2)0^(k3)05(k4) > 



It is seen that only matrix elements in Eqn. H29|) for which the indices a^-fS 
are pairwise equal survive disorder-averaging, and also that the average has 
no dependence on the energy of afi^S. In the spinless case, the first two terms 
on the right hand side make equal contributions and produce the constant 
charging energy in the Universal Hamiltonian of Eq. H26I) , while in the spinful 
case they produce the charging and exchange terms. The final term of Eq. 
H3UI) produces the Cooper interaction of Eq. (|26|) . 

Thus the UH contains the rotationally invariant part of the Landau in- 
teraction. The others, i.e., those that do not survive ensemble averaging, are 
dropped because they are of order 1/g. But we have seen before in the BCS 
instability of the Fermi liquid that a term that is nominally small to begin 
with can grow under the RG. That this is what happens in this case was 
shown by the RG calculation of Murthy and Mathur. There was however one 
catch. The neglected couplings could overturn the UH description for cou- 
plings that exceeded a critical value. However the critical value is of order 
unity and so one could not trust either the location or even the very existence 
of this critical point based on their perturbative one-loop calculation. Their 
work also gave no clue as to what lay on the other side of the critical point. 

Subsequently Murthy and I |24| showed that the methods of the large N 
theories (with g playing the role of N) were applicable here and could be 
used to show nonperturbatively in the interaction strength that the phase 
transition indeed exists. This approach also allowed us to study in detail the 
phase on the other side of the transition, as well as what is called the quantum 
critical region, to be described later. 



Oo,SO/3-,Okik4"k2k3 , Oo-y O^SOki k3 "k2 k4 
3 1 3 



'5Q/y<5-y5(5k^,-k2 (^ks .-k4 



(30) 
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Fig. 5. Feynman diagrams for the full four-point amplitude Fa/s-yS- 

Let us now return to Murthy and Mathur and ask how the RG flow is 
derived. After integrating some of the g states within Et, we end up with 
g' = ge~* states. Suppose we compute a scattering amplitude Fap-ys for the 
process in which two fermions originally in states aP are scattered into states 
7(5. This scattering can proceed directly through the vertex Vaf3-ys{t), or via 
intermediate virtual states higher order in the interactions, which can be 
classified by a set of Feynman diagrams, as shown in Figure All the states 
in the diagrams belong to the g' states kept. We demand that the entire 
amplitude be independent of g' , meaning that the physical amplitudes should 
be the same in the effective theory as in the original theory. This will lead 
to a set of flow equations for the Vafi-yS- In principle this flow equation will 
involve all powers of V but we will keep only quadratic terms (the one-loop 
approximation). Then the diagrams are limited to the ones shown in Figure 
leading to the following contributions to the scattering amplitude Fa^-ys 

Y^,i_NAt)_^V^p^.V,,,s (31) 

where the prime on the sum reminds us that only the g' remaining states are 
to be kept and where Npiot) is the Fermi occupation of the state a. We will 
confine ourselves to zero temperature where this number can only be zero 
or one. The matrix element Vap-fS now explicitly depends on the RG flow 
parameter t. 

Now we demand that upon integrating the two states at ±g' A/2 we re- 
cover the same Fap-^s. Clearly, since g' = ge~^, 
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The effect of this differentiation on the loop diagrams is to fix one of the 
internal lines of the loop to be at the cutoff ztg'A/2, while the other one 
ranges over all smaller values of energy. In the particle- hole diagram, since fj, 
or v can be at +g'A/2 or —g' A/2, and the resulting summations are the same 
in all four cases, we take a single contribution and multiply by a factor of 4. 
The same reasoning applies to the Cooper diagram. Let us define the energy 
cutoff A = g' A/2 to make the notation simpler. Since we are integrating out 
two states we have 5g' = 2 



dt 

fi—A.v '''' 

+ %4 E ''-^i-M;i^V^fs^.V,^^s (33) 

fi—A^iy 

where fi = A means e^, = A and so on. The changed sign in front of the 
1-loop diagrams reflects the sign of Eq. (^2)1 

So far we have not made any assumptions about the form of Vaf3-yS, and 
the formulation applies to any finite system. In a generic system such as an 
atom, the matrix elements depend very strongly on the state being integrated 
over, and the flow must be followed numerically for each different set afi'yS 
kept in the low-energy subspace. 

In our problem things have become so bad that are good once again: the 
wavefunctions (/!)(k) that enter the matrix elements above have so scrambled 
up by disorder that they can be handled by RMT. In particular it is possible 
to argue that the sum over the g' terms may be replaced by it ensemble 
average. In other words the flow equation is self-averaging. While the most 
convincing way to show this is to compute its variance, and see that it is 
of order l/^/g times its average, this fact can be motivated in the following 
way: There is a sum over 3> 1 values of with a slowly varying energy 
denominator, which makes the sum over v similar to a spectral average, 
which in RMT is the same as an average over the disorder ensemble. A more 
sophisticated argument is presented in Ref We can use the result 



< '/'^(Pi)C(P2)0>/(P3)?!'a'(P4) > = 
<5l4<523 S13S24 <^l,-2<53,-4 

g2 g3 g3 



(34) 



to deal with the product of four wavefunctions inside the loop and deal with 
the energy sum as follows: 



,A+\e,,\ J AA + e A 
^ 
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We are exploiting the fact that wavefunction correlations are energy inde- 
pendent in the large -g limit. 

After we make this simplifications we find that there are many kinds of 
terms of which one kind dominates in the large -g limit. 

Let us go back to the properly antisymmetrized matrix element defined 
in terms of the Fermi liquid interaction function, Eq. H29|l . Since there is a 
product of two y's in each loop diagram, and each V contains 4 terms, it is 
clear that each loop contribution has 16 terms. Let us first consider a term 
of the leading type in the particle-hole diagram, which survives in the large-g 
limit. Putting in the full wavefunction dependencies (and ignoring factors 
other than g, g') we have the following contribution from this type of term 



dt Leadmg ^'''"'kk'pp' 



x</.;(k)<^;(k')0^(k')05(k)0*(p)C(p')0.(p)04p')(36) 

Substituting the appropriate momentum labels for the particlc-hole dia- 
gram in Eqn. (|34|l . we see that the wavefunction average relevant to the sum 
over intermediate states is 

% - i±%P: (37) 
9 9 

Using the self-averaging, the first term of Ea. (|37|l forces p = p' in Eq. 
For large g, using 

p 

we obtain a convolution of the two Fermi liquid functions 

J2 - - ^k') = g{ul + 2 E "™ '^o^™^^ - ^')) (39) 

p m— 1 

where we have reverted to the notation = 9]^^, 9' = 9]^i. In the second term 
of Eq. H37|l . the (5p.-p' turns out to be subleading, while the other allows 
independent sums over p, p'. This means that only uo contributes to this 
term, (other avrerage to zero upon summing over all angles) which produces 

- J2 "(^k - 9p)u{9p, - 9^,) = g^ul (40) 
pp' 

Feeding this into full expression for this contribution to the particle-hole 
diagram, we find it to be 

= sLA\n2j: ( £ 
dt Leadmg ^ kk' V m=i 

C(k)0^(k')0^(k')05(k) (41) 
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Notice that the resuh is stiU of the Fermi hquid form. In other words the 
coupHngs Vapjs which were written in terms of Landau parameters m,„, flow 
into renormahzed couphng once again expressible in terms of renormahzed 
Landau parameters. By comparing the two sides, we see each Um flows inde- 
pendently of the others as per 

^ = -e-^iln2)ul m^O (42) 

The above equation can be written in a more physically transparent form 
by using a rescaled variable (for m ^ only) 

Urn = e^^Um (43) 

in terms of which the flow equation becomes 
du 

= -u,n - (In 2)ul = (3{u,n) (44) 

where the last is a definition of the /3-function. 

The reason Uo does not flow is that the corresponding interaction com- 
mutes with the one-body "kinetic" part, and therefore does not suffer quan- 
tum fluctuations. 

This is the answer at large g. We have dropped subleading contributions 
of the following type: 



dVap^s I 



v=-A 



kk' pp' 



xC(p)0^(k')0^(k')</>5(k)</);(k)C(p')0.(p)0p(p') (45) 



Note that the momentum labels of 0* and 0* have been exchanged com- 
pared to Eq. and there is a minus sign, both coming from the antisym- 
metrization of Eq. H29(l . Once again we ensemble average the internal /z, v 
sum, the wavcfunction part of which gives 



< <^;(k)<^:(p')0.(p)0,(p') ^-^^ 4p+4 p-^p,-p- 

g 5 

It is clear that there is an extra momentum restriction in each term com- 
pared to Eq. (|37|l . which means that one can no longer sum freely over p to 
get the factor of g in Eq. (|39|l , or the factor of g^ in Eq. H4U|) . Therefore this 
contribution will be down by 1 / g compared to that of Eq. H36|) . 

Turning now to the Cooper diagrams, the internal lines are once again 
forced to have the same momentum labels as the external lines by the Fermi 
liquid vertex, therefore they do not make any leading contributions. 
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The general rule is that whenever a momentum label corresponding to an 
internal line in the diagram (here fi and v) is forced to become equal to a 
momentum label corresponding to an external disorder label (here a,(3,"f, or 
5), the diagram is down by 1/g, exactly as in the 1/N expansion. The fact 
that 1/g plays the role of 1/iV was first realized by Murthy and Shankar 
j24| . Not only did this mean that the one loop flow of Murthy and Mathur 
was exact, it meant the disorder-interaction problem of the chaotic dot could 
be solved exactly in the large g limit. It is the only known case where the 
problem of disorder and interactions |26[ 1771 128| can be handled exactly. 

From Eqn. (|44|l it can be seen that positive initial values of Um (which 
are equal to initial values of u„i inherited from the bulk) are driven to the 
fixed point at Um = 0, as arc negative initial values as long as Um,(t = 0) > 

= — 1/ In 2. Thus, the Fermi liquid parameters are irrelevant for this range 
of starting values. Recall that setting all u,n = for m 7^ results in the 
Universal Hamiltonian. Thus, the range Um > ujj^ is the basin of attraction 
of the Universal Hamiltonian. On the other hand, for Um{t = 0) < results 
in a runaway flow towards large negative values of Um, signalling a transition 
to a phase not perturbatively connected to the Universal Hamiltonian. 

Since we have a large N theory here (with N = g), the one- loop flow and 
the new fixed point at strong coupling are parts of the final theory. ^What is 
the nature of the state for Um{t = 0) < m^? The formalism and techniques 
needed to describe that are beyond what was developed in these lectures, 
which has focused on the RG. Suffice it to say that it is possible to write 
the partition function in terms of a new collective field <t (which depends 
on all the particles) and that the action S{a) has a factor g in front of it, 
allowing us to evaluate the integral by saddle point (in the limit g 00), 
to confidently predict the strong coupling phase and many of its properties. 
Our expectations based on the large g analysis have been amply confirmed 
by detailed numerical studies |2SI- For now I will briefiy describe the new 
phenomenon in qualitative terms for readers not accustomed to these ideas 
and give some references for those who are. 

In the strong coupling region <t acquires an expectation value in the 
ground state. The dynamics of the fermions is affected by this variable in 
many ways: quasi-particle widths become broad very quickly above the Fermi 
energy, the second difference A2 has occasionally very large values and can 

^ However the exact location of the critical point cannot be predicted, as pointed 
out to us by Professor Peet Brower. The reason is that the Landau couplings 
defined at a scale El much higher than Et (but much smaller than Ep) 
and their flow till we come down to Et, where our analysis begins, is not within 
the regime we can control. In other words we can locate u* in terms of what 
couplings we begin with, but these are the Landau parameters renormalized in 
a nonuniversal way as we come down from El to Et- 
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even be negative, ^ , and the system behaves hke one with broken time-reversal 
symmetry if m is odd. 

Long ago Pomeranchuk [221 found that if a Landau function of a pure 
system exceeded a certain value, the fermi surface underwent a shape trans- 
formation from a circle to an non-rotationally invariant form. Recently this 
transition has received a lot of attentionPHU^ The transition in question is 
a disordered version of the same. Details are given in Refs. (|211, |2S1)- 

Details aside, there is another very interesting point: even if the coupling 
does not take us over to the strong-coupling phase, we can see vestiges of the 
critical point and associated critical phenomena. This is a general feature 
of many quantum critical points^^, i.e., points like where as a variable 
in a hamiltonian is varied, the system enters a new phase (in contrast to 
transitions wherein temperature T is the control parameter). 

Figure shows what happens in a generic situation. On the x-axis a vari- 
able {um in our case ) along which the quantum phase transition occurs. 
Along y is measured a new variable, usually temperature T. Let us consider 
that case first. If we move from right to left at some value of T, we will first en- 
counter physics of the weak-coupling phase determined by the weak-coupling 
fixed point at the origin. Then we cross into the critical fan (delineated by the 
V-shaped dotted lines), where the physics controlled by the quantum criti- 
cal point. In other words we can tell there is a critical point on the x -axis 
without actually traversing it. As we move further to the left, we reach the 
strongly-coupled symmetry-broken phase, with a non-zero order parameter. 

It can be shown that in our problem, 1/g plays the role of T. One way to 
see is this that in any large N theory N stands in front of the action when 
written in terms of the collective variable. That is true in this case well for 
g. (Here g also enters at a subdominant level inside the action, which makes 
it hard to predict the exact shape of the critical fan. The bottom line is 
that we can see the critical point at finite 1/g. In addition one can also raise 
temperature or bias voltage to see the critical fan. 

Subsequent work has shown, in more familiar examples that Landau in- 
teractions, that the general picture depicted here is true in the large g limit: 
upon adding sufficiently strong interactions the Universal Hamiltonian gives 
way to other descriptions with broken symmetry]^ 12] • 

^ How can the cost of adding one particle be negative (after removing the charging 
energy)? The answer is that adding a new particle sometimes lowers the energy 
of the collective variable which has a life of its own. However, if we turn a blind 
eye to it and attribute all the energy to the single particle excitations, A2 can 
be negative. 

^ The only result that is not exact in the large g limit is the critical value since 
the input value of Um at Et is related in a non-universal way to the Landau 
parameter. In other words, the Landau coupling tt„i settles down to its fixed 
point value at an energy scale that is much larger than Et- We do not know how 
it flows as we reach energies inside Et wherein our RMT assumptions are valid. 
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Fig. 6. The generic phase diagram for a second-order quantum phase transition. 
The horizontal axis represents the coupling constant which can be tuned to take one 
across the transition. The vertical axis is usually the temperature in bulk quantum 
systems, but is l/g here, since in our system one of the roles played by g is that of 
the inverse temperature. 
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